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EFFICIENT  SECOND  HARMONIC  CONVERSION  OF  BROADBAND 
HIGH-PEAK-POWER  Nd:GLASS  LASER  RADUTION  USING 
LARGE-APERTURE  KDP  CRYSTALS  IN  QUADRATURE 


1 .  Introduction 

Large  Nd:glass  laser  systems  have  been  developed  as  tools  for  laiser  fusion 
research.  Recent  laser-matter  interaction  experiments  have  shown  that  a  laser 
wavelength  shorter  than  one  micron  is  required  for  efficient  target  coupling  and 
suppression  of  some  laser-plasma  instabilities  [1] - [5] .  For  directly  driven 
fuel  pellets,  one  also  requires  nearly  uniform  illumination  to  produce  the 
highly-symmetric  ablation  pressure  necessary  for  high  gain.  Estimates  for  this 
symmetry  suggest  that  it  must  be  better  than  one  or  two  percent  [6] ,  [7] .  These 
requirements  have  stimulated  the  development  of  techniques  to  improve  focal 
uniformity  [8] -[11]  and  efficiently  produce  shorter  wavelength  light  by  harmonic 
conversion  [12] - [18] . 

One  of  the  more  promising  beam  smoothing  techniques,  Induced  Spatial 
Incoherence  (ISI)  [8] ,  [10] ,  requires  laser  radiation  with  a  short  coherence  time 
(r  <  2  ps)  and  therefore  a  broad  bandwidth  (Av/c  >  15  cm~^) .  Bandwidths  of  up 
to  30  cm'l  are  readily  achievable  in  Nd: glass  laser  systems  in  the  spectral 
region  near  one  micron.  Evaluation  of  the  ISI  technique  at  shorter  wavelength 
using  a  Nd: glass  laser  system  requires  a  harmonic  conversion  process  that  also 
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maintains  broad  bandwidth.  Here,  we  report  the  production  of  broadband  527-nffl 
has  several  advantages  over  the  conventional  single  Type  II  crystal.  In  a 
single  Type  II  crystal,  birefringence  limits  the  conversion  efficiency  by 
creating  a  group  velocity  mismatch  between  the  two  orthogonal  components  of  the 
fundamental.  The  quadrature  configuration  compensates  for  thjis  mismatch  because 
the  transverse  projections  of  the  two  crystal’s  principle  axes  are  orthogonal, 
as  shown  in  Fig.  1.  Hence,  the  fundamental  component  that  is  polarized  along 
the  slow  axis  in  the  first  crystal  will  be  polarized  along  the  fast  axis  in  the 
second.  This  arrangement  also  retains  the  major  advantages  of  Type  II  operation 
(i.e.,  a  wide  input  intensity  range  over  which  conversion  is  high  [18]  and  a 
relatively  high  tolerance  to  small  angular  misalignments  and  beam  divergence) . 
Furthermore,  harmonic  back-conversion  is  minimized  because  the  second  harmonic 
light  produced  in  the  first  crystal  is  not  at  the  correct  polarization  for 
interaction  in  the  second. 

A  theoretical  model  which  simulates  the  broadband  conversion  process  is 
described  in  Section  II.  Using  this  model,  we  calculated  the  conversion 
efficiency  and  spectral  properties  of  the  harmonically  generated  radiation  for 
two  specific  input  bandwidths.  Details  of  the  experimental  arrangement  are 
presented  in  Section  IH.  The  experimental  results  are  presented  in  Section  IV 
where  reausonable  agreement  with  the  numerical  calculations  is  also  shown. 
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II.  Theory  of  Broadband  Eamonic  Conversion 

Broadband  harmonic  conversion  can  be  described  by  the  usual  theory  of  three 
wave  mixing  in  dispersive  birefringent  nonlinear  crystals  [20] .  Each  of  the 
colinearly  propagating  waves  {i(a,t),  can  be  written  in  the 

form 


{n(*»^)  =  ^/2En(*»t)exp[i(kn8  -  «nt)]  +  c.c.  ,  (1) 

where  the  spatial  and  temporal  angular  carrier  frequencies  are  chosen  to  satisfy 
the  exact  phase-matching  conditions 

ki  k2  -  k3  =  0  ,  11^1  +  *>2  -  ^3  =  0  ,  (2) 

and  the  complex  amplitudes  £Q(z,t)  are  assumed  to  be  slowly -varying  on  the  scale 
of  kjj'^  and  For  the  broadband  case,  these  amplitudes  must  account  for  the 

statistical  properties  of  the  light,  plus  any  phase  mismatch  effects  on  spectral 
components  that  are  detuned  from  the  carrier  frequencies. 

The  coupled  amplitude  equations  are  [21] , [22] 

(d/8z  +  ai8/8t  +  ibi82/8t2  +  ^i)Ei  =  -  W1CB3E2*  ,  (3a) 

(8/8z  +  a28/8t  +  ib282/8t2  +  /J2)E2  =  -  ^2^3®!*  » 
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(d/8s  +  ib382/0t2  +  Pz)^3  ~  W3CB1E2  ,  (3c) 

where  Eq  =  Vgn“l  -  Vg3“^  ,  Vgn”^  =  ^kn/Bw^  is  the  inverse  group  velocity  of  the 
nth  wave,  bjj  =  accounts  for  the  group  velocity  dispersion,  is 

the  linear  loss  term,  and  C  is  the  nonlinear  coupling  constant.  If  the 
amplitudes  are  scaled  as  IEqI^  i  Iq  ,  where  Iq  is  the  intensity  in  W/cm^  and  the 
lengths  and  velocities  in  equations  (3)  are  in  units  of  cm  and  cm/s 
respectively,  then  the  coupling  constant  is 

C  =  (2  X  106/eoc3nin2n3)l/2(d/Co)sin(mBa)  , 

where  Cq  is  the  permittivity  of  free  space  in  MKS  units,  c  is  the  velocity  of 
light,  and  nQ  are  the  refractive  indices.  The  nonlinear  optical  coefficient  is 
defined  such  that  for  KDP  d/£Q  =  0.39  pm/Y.  The  phase  matching  angle  is  Bq,;  for 
Type  I  phase  matching  m  =  1,  and  for  Type  II  m  =  2.  Equations  (3)  describe  the 
interaction  as  seen  by  an  observer  moving  along  with  the  £3  amplitude.  For  most 
cases  of  interest,  the  contributions  due  to  group  velocity  dispersion  are  small 
in  comparison  to  the  other  terms. 

In  second  harmonic  conversion,  E^  *  Etfx  and  E2  *  E^/y  represent  orthogonal 
spatial  components  of  the  fundamental  amplitude  E^  =  ^2  ^  >  vliHe  E3  £2^ 

describes  the  harmonic  (W3  =  2u) .  The  limitations  on  2u  conversion  imposed  by 
finite  optical  bandwidth  can  be  easily  illustrated  in  the  limit  of  negligible 
pump  depletion  IE2ut  «  lE^I.  If  one  also  ignores  the  small  absorption  and 
group  velocity  dispersion  effects,  then  the  mixing  equations  reduce  to 
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(a/Bz  +  ai8/at)E,«  s  0  , 


(4a) 


(a/8z  +  a28/8t)E|^  -  0  , 


(4b) 


8E2(„/8z  2  2wCEy3jEyy 


(4c) 


For  a  crystal  of  thickness  L  and  the  usual  boundary  condition  ^2u  =  0  at  z  =  0, 
the  solution  is 


fL 

E2«(t)  =  2tfC  1  E^(t  -  aiz)Eh^y(t  -  a2z)dz 

Jo 


(5) 


It  is  instructive  to  first  consider  the  case  shere  the  fundanental  is  a 
Bonochromatic  wave  whose  actual  frequency  is  detuned  from  u  by  some  small  amount 
Su;  i.e.,  E|^ j  (t)  *♦  Ajexp(-i5wt) ,  where  Aj  ♦  Ajj,  Ay  are  constants.  Then 


E2tf('b)  *  2tfCAxAyLexp[-2iffw(t  -  l/2AVj.j.)]sinc(5w/AV{.r)  i 


(6) 


where  sinc(x)  =  sin(x)/x  and 


A*'cr  = 


2/L  1/L 

*1  +  a2  ^/2(vgr^  +  ^g2"^)  -  ▼gS’^ 


(7) 


is  the  frequency  bandpass  (Atfcr/Sjr)  of  the  crystal  in  the  low  conversion  limit 
[23]. 


5 


la  the  broadband  case,  the  incident  fundamental  is  modeled  by  quasi-stationary 
chaotic  noise  [24] .  Its  spectrum  is  comprised  of  many  randomly-phased  modes  of 
total  bandwidth  centered  around  w/2r.  Thus  £(,r(t)  is  a  stochastic  function 
characterized  by  a  coherence  time  t^  =  l/Ai/p  much  shorter  than  the  laser  pulse 
duration.  Applying  these  considerations  to  expression  (5) ,  .ope  observes  that  if 
la^lL,  la2lL  «  t(.  ,  then  E2u(t]  is  given  by  the  phase-matched  result 
2tfCE(^(t)E^(t)L.  However,  if  either  laxIL  or  la2lL  become  comparable  to  tc 
the  integrand  can  undergo  a  phase  shift  within  the  (0,L)  interval,  resulting  (on 
the  average)  in  lower  harmonic  conversion.  The  interpretation  of  these  results 
is  straightforward  in  Type  I  operation,  where  ax  =  a2  -  aj.  According  to 
expression  (7),  the  criterion  lajlL  «  t(.  then  reduces  to 

Ai/p  «  Ai/cr  ;  (8) 

i . e . ,  the  optical  bandwidth  must  be  much  less  than  the  crystal  bandpass .  For 
1.054  pm  light  in  KDP,  we  find  aj  =  0.0529  ps/cm,  thus  giving  -  630  cm'^ 

(280  cm-1  FTiraM)  for  L  a  1  cm.  Type  I  operation  offers  the  advantages  of 
simplicity  and  broad  spectral  bandpass,  but  the  angular  tolerance  and  the 
intensity  range  over  which  the  conversion  remains  high  both  tend  to  be 
relatively  narrow  in  comparison  to  those  of  Type  II  operation. 

In  Type  II  operation,  the  conversion  may  be  limited  more  by  the  crystal’s 
birefringence  than  by  its  bandpass.  Consider  for  example  Type  II  operation  of 
KDP  at  1.054  pm.  If  E^x  lies  along  the  o  axis  and  E^^y  lies  in  the  oe  plane, 
then  ax  **  ag  =  0.5539  ps/cm  and  a2  **  a^  =  -0.7669  ps/cm.  Expression  (7)  then 
gives  Al/gr/c  =  313  cm~^  (138  cm~l  FWEH)  for  L  =  1  cm,  which  suggests  that  high 
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conrersion  could  be  achieved  ia  a  single  crystal  with  the  20  ca~^  to  30  cb~^ 
optical  bandeidtbs  available  from  Ndtglass  lasers.  However,  criterion  (8)  is  no 
longer  sufficient  to  ensure  high  conversion  because  aj^  and  a2  sre  of  opposite 
sign.  An  observer  moving  along  z  with  the  £21^  amplitude  (as  described  in 

(^))  would  see  £^  retarded  by  ag  ps/cm,  while  £|^  advanced  by  la^l  ps/cm. 
These  amplitudes  can  remain  correlated,  thus  allowing  the  integrand  of  £q.  (5) 
to  remain  constant  along  z,  only  if  the  condition 

hcA  >  *0  *  =  1.321  ps/cm  (9) 

is  satisfied.  In  a  single  crystal,  this  criterion  will  effectively  limit  the 
intensity  conversion  to  ~  50%  at  bandwidths  of  20  cm~^  to  30  cm~^. 

Higher  conversions  can  be  achieved  without  sacrificing  the  advantages  of  Type 
n  operation  by  using  the  quadrature  configuration  [18]  shown  in  Fig.  1.  The 
first  crystal  is  oriented  as  described  above.  At  the  output,  £^  is  retarded  by 
SgL  while  E|^  is  advanced  by  laglL.  This  crystal  can  convert  30%  to  50!S  of  the 
light  at  bandwidths  of  20  cm*^  30  cm~^  if  (ag  lagOL  **  tg.  In  the  second 
crystal  the  optic  axis  is  rotated  90*  around  the  propagation  direction,  thus 
interchanging  the  transverse  projections  of  the  o  and  e  axes.  The  amplitude  Ewx 
will  now  begin  to  advance  at  the  rate  la^l  ps/cm,  while  retards  at  the  rate 
ag  ps/cm.  If  this  crystal  has  a  thickness  L'  =  2L,  then  Zyx  und  Ej,^  will  be 
back  in  step  around  its  midplane,  and  its  harmonic  contribution  £2^’  can  be 
generated  efficiently  along  most  of  its  path.  The  harmonic  contribution  £2^ 
generated  by  the  first  crystal  is  orthogonal  to  E2v’-  Thus  £21,^  is  not  phase- 
matched  in  the  second  crystal,  and  will  propagate  through  it  with  only  the 
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lioear  losses  x£2tf  due  to  absorption  and  Fresnel  reflection.  The  net  intensity 
conversion  is  then  given  by 


(1  -  k'^)<12u>  *  <l2w*> 

,  - -  ,  (10) 

<I«(0,t)> 

where  <In>  =  <IEn|2>  is  the  intensity  of  the  nth  wave  averaged  over  times  much 
longer  than  t^  and  <1^(0, t)>  is  the  total  fundamental  <'Lux>  *  <Xuy>  incident  at 
the  first  crystal.  Using  low  loss  crystals  and  high  quality  AR  coatings,  one 
can  achieve  Intensity  conversions  greater  than  75%  with  output  bandwidths  of 
•*  30  cm~^. 

Broadband  2ii>  conversion  has  been  modeled  in  detail  by  use  of  a  numerical  code 
FAST2F  to  integrate  Eqs.  (3).  The  temporal  variation  is  handled  by  fast  Fourier 
transform  techniques  similar  to  those  described  elsewhere  [21] , [25] .  If  En(z,fl) 
is  the  Fourier  transform  of  En(*,t),  i.e. 

En(*»D)  =  TF{En(s,t)}  ,  Ea(£,t)  =  Trl{En(*,Q)}  ,  (11) 

then  Eqs.  (3)  take  the  general  form 

(a/5z  -  ianO  -  ibnflZ  ^  ,  (12) 

where  Mq  is  the  mixing  term  for  the  nth  wave.  (For  example,  Ui  =  -  W1CE3E2*.) 
The  group  velocity  tenss,  which  could  cause  a  rapid  variation  at  larger  values 
of  0,  are  removed  from  the  integration  by  using  the  unitary  transformations 
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En(*.Q)  =  En'expCiaafl  +  ibnQ2) 


(13) 


to  rewrite  £q.  (12)  in  the  'interaction  representation,' 

(B/8z  +  /yjEn*  =  exp(-ianQ  -  ibnQ2)Tp{Mn},  ,  (14) 

Equations  (14)  are  then  integrated  by  a  predictor-corrector  technique,  using 
Eqs.  (11)  and  (13)  at  each  new  point  along  z.  This  algorithm  was  benchmarked 
against  the  well-known  analytic  result  [26]  for  monochromatic  light  detuned  from 
the  exact  phase  matching  frequency;  the  agreement  was  better  than  1%  over  a  wide 
range  of  conversion  efficiencies. 

FAST2F  calculates  the  total  harmonic  conversion  and  the  spectral  and 
statistical  properties  of  the  2w  light;  it  can  be  configured  for  either  Type  I 
or  quadrature  Type  II  operation.  It  also  evaluates  the  correlation  between  the 
harmonic  amplitudes  generated  by  two  correlated  but  physically  separate  input 
beams  differing  only  in  their  intensities.  With  that  data,  one  can  assess  the 
possible  deterioration  of  transverse  coherence  (and  hence  bean  quality)  that  can 
sometimes  occur  when  a  crystal  is  strongly  driven  by  a  spatially  nonuniform 
input  intensity  [27].  For  the  simulations  reported  here,  these  hairmonic 
amplitudes  remained  correlated  to  better  than  60%.  The  crystal  parameters  used 
for  the  Type  II  simulations  are  listed  in  Table  I.  All  of  the  refractive  index 
and  group  velocity  parameters  are  taken  from  data  of  Zernike  [28] ;  the 
reflectivities  of  the  A£  coatings  (described  in  the  next  section)  and  the 
internal  loss  coefficients  were  specified  by  their  respective  manufacturers 
[29]. 
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The  iocident  ftmdaaental  was  aodeled  bj  an  array  of  Gaussian-distributed 
randoB  complex  numbers  whose  power  spectrum  envelope  matched  the  time-integrated 
laser  spectra  measured  in  the  experiment  (see  Section  lY) .  This  resulted  in 
stochastic  temporal  behavior  characterised  by  coherence  time 

tc  s  1/Ai/p,  which  remained  statistically-stationary  over  the  {entire  calculation 
range  (0,T)  »  t(..  In  order  to  avoid  aliasing  behavior,  the  intensity  was 
terminated  within  a  short  transition  interval  t^^  (tQ  <  t^j.  «  T)  at  both  ends 
of  the  calculation  range.  For  the  simulations  reported  here,  the  width  of  the 
remaining  (statistically-stationary)  interval  was  Tg  =  580  ps,  while  the 
coherence  times  were  typically  "  1  ps.  All  intensity  averages,  such  as  those 
used  in  £q.  (10),  were  evaluated  over  this  interval. 

Fig.  2  shows  a  typical  simulation  of  the  stochastic  temporal  behavior,  power 
spectrum,  and  autocorrelation  function  for  the  incident  fundamental  where  the 
intensity  averaged  over  the  sampling  interval  is  1  GW/cm^.  In  the  spectral 
plot,  the  power  is  averaged  over  *  1  cm*^  in  order  to  roughly  model  the  limited 
resolution  of  the  spectrograph  used  to  make  the  measurements.  The  2u 
simulations  were  carried  out  assiuning  that  the  crystals  are  phase-matched  to  the 
center  of  the  incident  fundamental;  hence  the  zero  frequency  was  chosen  to  lie 
at  the  centroid  of  this  spectr\un.  The  estimated  FWHM  width  of  the  spectrum  is 
17  CB~^.  From  the  autocorrelation  plot,  we  calculate  a  coherence  time 
t^  -  2  ps;  the  corresponding  spectral  width  l/t(.  *  16.7  cm~^  is  consistent  with 
the  estimated  FWHU. 

Figs.  3  and  4  show  the  calculated  temporal  behavior,  spectra,  and 
autocorrelation  functions  for  the  2u  generated  by  each  of  the  quadrature 
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crystals.  lo  this  example,  the  first  crystal  converted  ~  43%  of  the  incident 
fundamental  (described  in  Fig.  2) ,  while  the  second  converted  an  additional 
29%.  The  FVHM  bandwidths  estimated  from  the  spectra  are  again  consistent  with 
the  effective  values  l/t^  found  from  the  autocorrelation  functions,  (i.e., 

26  cm”^  FWBU  vs.  l/t,.  28  cm“^  for  the  first  crystal,  and  SS,  cm~^  FWEK  vs. 

1/tc  33  cm~^  for  the  second.)  The  bandwidths  from  the  first  crystal  generally 
remained  within  1.4  to  1.6  times  that  of  the  fundamental,  but  the  spectra  from 
the  second  crystal  broadened  significantly  at  higher  intensities,  as  shown  in 
Fig.  5.  This  appears  to  result  from  a  broadening  of  the  fundamental  by  the 
conversion  process  in  the  first  crystal;  at  higher  intensities  that  process  will 
preferentially  deplete  the  central  portion  of  the  fundamental  spectrum.  In 
practice,  this  effect  will  result  in  only  a  modest  broadening  of  the  total 
harmonic  spectrum  because  most  of  the  conversion  at  high  intensities  occurs  in 
the  first  crystal.  The  harmonic  components  generated  by  each  crystal  are 
orthogonally  polarized,  so  the  total  second  harmonic  spectrum  may  be  calculated 
by  simply  adding  the  individual  spectral  intensities.  Fig.  6  shows  the 
calculated  total  harmonic  spectrum  for  an  incident  fundamental  with  1  GW/cm^ 
average  intensity  and  17  cm~^  bandwidth.  At  this  intensity,  where  the  2u 
spectrum  is  dominated  by  high  conversion  in  the  first  crystal,  the  estimated 
FWBU  bandwidth  is  28  cm'^. 


Fig.  7  shows  calculated  intensity  conversion  efficiencies  for  two  incident 
fundamental  bandwidths .  The  results  were  evaluated  using  expression  (10) ,  with 
all  intensities  averaged  over  the  580  ps  interval  discussed  above.  In  order  to 
estimate  the  accuracy  of  the  results,  we  ran  the  simulations  for  the  0.5  GW/cm^ 
and  2.0  GW/cm^  points  using  four  statistically-independent  realizations  of  the 
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random  fundamental  amplitude  at  the  input.  The  RkS  deviation  from  the  average 
values  was  ~  0.5%  for  the  total  conversion  (solid  curves)  and  **  0.7%  for  the 
first  crystal  (dashed  curves) .  It  is  apparent  from  these  curves  that  quadrature 
conversions  are  not  only  higher  than  those  of  a  single  crystal,  but  they  also 
remain  high  over  a  much  wider  range  of  input  intensities.  In^  fact,  this  feature 
is  one  of  the  most  important  advantages  of  the  quadrature  configuration  [18] . 

The  curves  also  show  that  there  is  little  penalty  for  using  the  broader 
bandwidth  (27  cm~^)  fundamental  to  achieve  output  bandwidths  in  excess  of 
35  cm~^.  For  the  2.5  ns  (FWHU)  pulses  used  in  the  experiment,  the  peak 
intensity  was  limited  to  ~  1  GW/cm^  by  the  ~  5  J/cm^  damage  threshold  of  the 
crystals.  The  curves  indicate  that  higher  conversions  should  be  possible  with 
the  high  damage  threshold  crystals  currently  under  development  at  Osaka 
University  [30] . 


III.  Experimental  Apparatus 

The  fundamental  light  for  the  conversion  efficiency  measurements  was  produced 
by  the  PHAROS  IH  laser  system.  The  laser  (see  Fig.  8)  is  an  image-relayed,  Nd- 
doped  phosphate  glass  system  developed  at  NRL  to  study  laser-matter  interaction. 
The  laser  produces  up  to  ~  800  J  per  beam  of  1.054  ftm  light  at  the  output  of  the 
final  15-cm  clear  aperture  disc  amplifier.  The  beam  is  expanded  by  a  telescope 
to  a  20-cm  diameter  before  frequency  conversion.  For  this  investigation,  the 
maximum  energy  incident  on  the  crystals  was  **  500  J  in  a  2.5  ns  pulse  with 
typical  FWHM  bandwidths  of  either  ~  17  cm~^  or  **  27  cm~^.  The  corresponding 
maximum  peak  intensity  averaged  across  the  beam  was  **  1  GW/cm^. 
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The  bandwidth  of  the  laser  system  was  varied  by  changing  the  bandwidth  of  the 
oscillator.  The  broad  fluorescence  linewidth  of  the  Nd-phosphate  glass 
amplifier  used  in  the  oscillator  allows  many  longitudinal  modes  to  exist  in  the 
60-cm  long  cavity.  Our  oscillator  (shown  in  Fig.  6)  produces,  a  dO-ns  long  pulse 
(FWHM)  in  the  TEUqq  mode.  The  oscillator  can  be  operated  with  or  without  the 
angle-tuned  etalon  shown  in  the  figure.  A  typical  output  spectrum  with  an 
estimated  FWHM  of  15  cm~^  is  shown  in  Fig.  10a  where  the  oscillator  was  operated 
without  the  intra-cavity  etalon.  The  structure  in  the  spectrum  is  attributed  to 
weak  etalon-like  behavior  of  the  stacked-plate  polariser.  Broader  bandwidths 
were  obtained  using  a  0.10-mm  thick  etalon  with  a  reflectivity  of  1.7!(  per 
surface  oriented  in  the  cavity  to  suppress  the  center  of  the  output  spectrum 
producing  the  lineshape  shown  in  Fig.  10b.  A  wide  variety  of  spectral 
lineshapes  can  in  fact  be  produced  using  various  combinations  of  etalons  and  by 
controlling  the  round-trip  cavity  gain. 

A  nominal  2.S-n8  long  Gaussian-like  pulse  was  sliced  from  the  peak  of  the 
longer  oscillator  pulse  using  two  Fockel’s  cells  in  series  driven  by  a  fast 
("’  1  ns  risetime) ,  high-voltage  Kry tron  pulse  generator .  The  contract  ratio  of 
the  Pockel’s  cell  pulse  slicer  is  greater  than  10^.  The  temporally  shaped  pulse 
is  then  amplified  before  being  split  into  three  equal  intensity  beams  for 
delivery  to  the  large-aperture  amplifier  chains.  Each  amplifier  chain  consists 
of  a  double-passed  45-mm  diameter  rod  amplifier  and  three  disc  amplifiers 
ranging  in  clear  aperture  from  66  mm  to  150  mm.  Image-relaying  telescopes 
transport  and  expand  the  beam  from  an  initial  45-mm  diameter  to  the  final  150-mm 
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diameter.  In  addition/  Pockel’s  cells  and  Faradax  Botators  are  used  for 
interstage  isolation  and  protection  from  target  backscatter. 

The  degree  of  polarisation,  spatial  quality  and  divergence  of  the  fundamental 
beam  can  significantly  effect  the  harmonic  conversion  efficiency.  A  major 
concern  in  high-peak- poser  glass  laser  systems  is  the  effect  on  beam  quality  of 
high-intensity  spatial  fluctuations  enhanced  by  nonlinear  phase  shifts.  In 
addition  to  damaging  optics,  these  effects  may  also  cause  intensity  dependent 
depolarisation  and  beam  break-up,  as  has  been  observed  elsewhere  [31] , [32] . 
Image-relaying  and  spatial  filtering  can  ameliorate  but  not  eliminate  these 
problems.  PHAROS  III  is  image-relayed  but  not  tightly  spatially  filtered;  thus 
there  was  some  concern  that  beam  quality  might  limit  the  harmonic  conversion 
efficiency.  An  attempt  to  measure  the  depolarisation  of  the  beam  was  made  using 
a  large-aperture  thin-film  polariser  and  calorimeters  to  monitor  the  incident 
and  rejected  energy.  The  results  sere  limited  by  the  resolution  and  accuracy  of 
the  calorimeters,  but  they  indicate  that  the  beam  is  linearly  polarised  to 
better  than  a  fes  percent.  The  spatial  quality  of  the  beam  sas  assessed  by 
imaging  an  attenuated  reflection  of  the  beam  onto  Polaroid  positive/negative 
film.  To  calibrate  the  film,  the  beam  sas  passed  through  a  sedged  rattle  plate 
to  produce  multiple  images  on  the  film  which  differed  in  intensity  by  a  known 
amount.  The  negatives  sere  scanned  through  the  center  line  with  a  densitometer 
to  produce  film  density  curves.  The  resulting  curves  are  nearly  linearly 
proportional  to  the  time-integrated  intensity.  A  typical  profile  is  shown  in 
Fig.  11.  The  densitometer  scan  taken  along  one  diameter  shows  a  nominally  flat- 
topped  intensity  distribution  with  local  spatial  nonuniformities  of  less  than 

i 
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*  sox  about  tbe  average.  Tbe  divergence  of  the  beam  was  inferred  from  the 
diffraction  pattern  produced  by  passing  tbe  beam  through  a  mask  containing  an 
array  of  circular  apertures.  Using  this  method,  we  measured  the  maximum  full- 
angle  divergence  of  the  beam  to  be  130  /(rad.  Measurements  of  the  focal 
diameter  indicated  that  the  beam  was  **  7  times  tbe  diffractiop  limit.  We 
concluded  from  these  results  that  the  beam  was  of  adequate  quality  to  perform 
the  conversion  efficiency  measurements  without  any  additional  effort  to  improve 
it. 

The  quadrature  harmonic  conversion  assembly  is  a  monolithic  structure 
consisting  of  an  aluminum  frame  in  which  two  25-cm  square  KDP  crystals  are 
mounted.  The  thicknesses  of  the  first  and  second  crystals  are  1.0  cm  and 
2.5  cm,  respectively.  Both  of  the  crystals  are  dual-wavelength,  anti-reflection 
coated  on  both  sides  using  a  Sol-Gel  process  [33] .  The  surface  reflectivity 
using  this  process  is  typically  »  2X  per  surface  at  lu  and  ~  IX  per  surface  at 
2u  [34] .  Direct  application  of  the  coating  to  the  crystal  surface  alleviates 
the  Raman  scattering  problem  that  can  occur  when  index  matched  fluid  cells  are 
used  [15] .  The  crystals  and  their  mounting  frame  were  precision  machined  so 
that  the  optimum  propagation  axes  through  the  crystals  are  within  50  /(rad  of 
each  other  when  mounted  in  the  holder. 


IV.  Experimental  Results 

The  measurements  were  made  with  the  diagnostic  arrangement  shown  in  Fig.  12. 
Volume  absorbing  calorimeters  were  used  to  measure  tbe  energy  at  both  the 
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fundamental  and  second  harmonic  wavelengths.  Because  the  energy  measurement  is 
tine  integrated  and  harmonic  conversion  is  intensity  dependent,  both  the 
temporal  pulse  shape  and  the  spatial  profile  of  the  fundamental  are  required  to 
evaluate  the  energy  conversion  efficiency.  The  temporal  profile  was  measured 
with  a  silicon  FIN  diode  which  has  an  impulse  response  of  ~  .300  ps.  Ueasurement 
of  the  spatial  profile  was  described  in  the  previous  section.  Within  the 
resolution  of  the  measurements,  neither  the  temporal  pulse  shape  nor  the  spatial 
profile  varied  appreciably  from  shot  to  shot.  Thus  to  evaluate  the  conversion 
efficiency,  a  generic  temporal  pulse  shape  and  spatial  profile  were  convolved 
with  the  appropriate  theoretical  data  presented  in  Fig.  7.  This  results  in 
energy  conversion  efficiency  data  expressed  as  a  function  of  the  incident  energy 
for  a  specific  temporal  pulse  shape  and  spatial  profile.  In  Fig.  13,  this  data 
is  compared  to  energy  conversion  measurements  for  both  fundamental  baindwidth 
cases.  These  results  describe  the  external  conversion  efficiency  of  the  crystal 
pair.  (i*e.,  no  corrections  were  made  for  crystal  absorption  or  Fresnel  losses 
in  the  experimental  data,  but  were  included  in  the  theory.)  From  the  figure, 
one  can  see  good  agreement  between  theory  and  experiment.  Energy  conversion 
efficiencies  of  ~  55%  were  observed  for  laser  pulse  energies  of  ~  500  J  for  the 
17  cm~^  bandwidth.  The  conversion  efficiency  decreased  by  15%  (from  ~  55%  to 
~  47%  at  500  J)  when  the  bandwidth  was  increased  by  ~  60%.  For  comparison,  a 
15%  to  20%  increase  in  the  efficiency  (relative  to  the  17  cm~^  bandwidth  case) 
was  observed  when  a  time-bandwidth  limited  oscillator  was  used  with  the  laser. 

The  time-integrated  spectra  were  measured  with  two  0.25-a  spectrographs  having 
resolutions  of  0.05  nm  and  0.025  nm  for  the  fundamental  and  second  harmonic, 
respectively.  Typical  spectral  data  are  presented  in  Figs.  14  and  15. 
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Propagation  of  broadband  radiation  through  a  glass  anplifier  is  complicated  by 
the  limited  bandpass  of  the  glass  and  self-phase  modulation  due  to  the  nonlinear 
index.  The  former  tends  to  narrow  the  spectral  width  while  the  latter  broadens 
it.  For  our  operating  conditions,  it  is  not  known  whether  (or  to  what  degree) 
these  phenomena  are  contributing  to  spectral  modification;  however,  comparison 
of  the  fundamental  spectra  to  the  corresponding  oscillator  spectra  presented  in 
Fig.  10  shows  that  these  operating  conditions  allowed  propagation  through  the 
glass  amplifiers  with  little  effect  on  the  spectral  shape.  The  second  harmonic 
spectra  in  each  case  show  an  increase  in  bandwidth  by  a  factor  of  1.5  to  1.6 
compared  to  the  fundamental  (consistent  with  the  calculations  presented  in 
Section  II).  For  the  narrower  bandwidth  case,  one  finds  the  estimated 
experimental  width  of  **  2S  cm~^  for  the  second  harmonic  to  agree  well  with  the 
calculated  value  of  ~  28  cm~^;  differences  in  fine  structure  are  due  to 
different  statistical  realizations,  averaging  times  (2.5  ns  vs.  the  theoretical 
value  of  580  ps)  and  limited  resolution  in  the  measurement. 


Y.  Conclusions 


We  have  shown  that  second  harmonic  light  can  be  efficiently  produced  with 
sufficient  bandwidth  for  ISI  despite  the  inherent  narrowband  resonance  of  the 
phase-matched  harmonic  crystals.  Conversion  efficiencies  of  up  to  S5%  were 
achieved  for  fundamental  bandwidths  of  ~  17  cm~^  where  the  incident  intensities 
were  limited  by  the  long  laser  pulses.  At  the  broauler  27  cm~^  bandwidth, 
conversion  efficiency  was  only  moderately  reduced.  The  results  agree  well  with 


17 


code  calculations,  which  indicate  that  higher  conversion  efficiencies  are 
possible  by  increasing  the  fundamental  input  intensity. 


One  disadvantage  to  the  quadrature  conversion  scheme  is  the  random 
polarisation  of  the  harmonic  beam.  Some  plasma  instabilities,  are  sensitive  to 
polarisation;  thus,  a  defined  polarisation  on  target  is  preferred.  In 
principle,  it  is  possible  to  detune  the  first  crystal  in  the  quadrature 
configuration  while  maintaining  adequate  alignment  on  the  second,  even  though 
the  crystals  are  rigidly  mounted  together.  In  this  situation,  the  second 
crystal  becomes  the  primary  harmonic  converter  while  the  first  acts  as  a  wave 
plate,  thus  maintaining  the  advantages  of  the  quadrature  configuration. 
Preliminary  experimental  results  have  shown  that  when  the  first  crystal  is 
detuned  sufficiently  to  reduce  its  conversion  efficiency  by  a  factor  of  about 
fifteen  at  low  power,  the  resulting  conversion  efficiency  at  high  power  is 
reduced  only  by  *•  1055  for  the  17-cm~l  bandwidth  case. 
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Fig.  2.  Sioulations  of  the  incident  fundafflental  light  (1.054  urn 
the  stochastic  temporal  behavior  (average  intensity  is  1  GW/cm^) 
spectrum,  and  (c)  the  squared  autocorrelation  function.  , 


) ,  showing  (a) 

,  (b)  the  power 
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Fig.  3.  Simulations  of  the  second  harmonic  light  generated  in  the  first 
crystal,  showing  (a)  the  stochastic  temporal  behavior,  (b)  the  power  spectrum, 
and  (c)  the  squared  autocorrelation  function. 
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Fig.  5.  Calculated  bandwidths  of  the  harmonic  light  generated  in  the  first  and 
second  crystals  (dashed  and  solid  lines,  respectively)  vs.  the  average  intensity 
of  the  incident  fundamental  at  two  different  bandwidths. 
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Fig.  6.  Calculated  spectrum  of  the  total  second  harmonic  light  generated  by  tw 
crystals  in  quadrature  for  an  incident  fundamental  with  an  average  intensity  of 
1  GW/cm^  and  a  bandwidth  of  17  cm~^. 
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FREQUENCY  | - 1 

(W  9.2  cm-'' 

Fig.  10.  Typical  tine-integrated  oscillator  spectra  (center  wavelength  ~  1.054 
foa)  for  oscillator  operation  (a)  without  an  intra-cavity  etalon  and  (b)  with  an 
etalon. 
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Fig.  11.  Typical  spatial  profile  and  densitometer  scaui. 
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Fig.  12.  Diagnostic  arrangement  used  to  measure  the  second  harmonic  conversion 


(%)  AON3IOIdJ3  NOISa3ANOO 
OINOWUVH  aN003S 


35 


Fig.  13.  Second  harfflonic  energy  conrersion  efficiency  measured  for  the 
quadrature  crystal  arrangement.  The  theoretical  predictions  for  fundamental 
FWHk  bandwidths  of  17  cm“^  and  27  cm“^  are  shown  as  the  solid  and  dashed  lines, 
respectively.  The  peak  intensity  averaged  across  the  beam  was  ~  1  GW/cm^  for  a 
laser  pulse  energy  of  500  J. 
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Fig.  14.  Typical  spectra  for  (a)  the  fundanental  and  (b)  the  second  harnonic 
when  no  etalon  is  used  in  the  oscillator. 
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Fig.  15.  Typical  spectra  for  (a)  the  fundasental  and  (b)  the  second  harmonic 
when  an  etalon  is  used  in  the  oscillator  to  broaden  the  bandwidth. 
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